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Abstract 

This paper is devoted to a description of integrable Hamiltonian hydrodynamic 
chains associated with Dorfman Poisson brackets. Three main classes of these hy- 
drodynamic chains are selected. Generating functions of conservation laws and 
commuting flows are found. Hierarchies of these Hamiltonian hydrodynamic chains 
are extended on negative moments and negative time variables. Corresponding 
three dimensional quasilinear equations of the second order are presented. 
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1 Introduction 

Recently (see [24]), integrable hydrodynamic chains written in the conservative form 

dth k = d x fk(ho,hi,...,hk+i), k = 0,1,2,... (1) 

were completely described (this classification problem was established in [25], the integra- 
bility criterion based on the differential geometric approach was suggested in [9], where 
first two expressions fo{h ,hi) and f 2 {h Q , hi,h 2 ) were determined). Also, it was proved 
that these conservative hydrodynamic chains can be written via special coordinates, i.e. 
the so called moments A k (h , hi, hk) such that 

A k = fiA k x +1 + f A k x + A k+1 (s A° x + Sl Al) + A k {r A° x + nA x ) (2) 

+k[A k+1 {w A° x + Wl Al + w 2 A 2 x ) + A k (v A° x + v x A x x + v 2 A 2 x ) + A k -\u A x + Ul Al + u 2 A 2 x )], 

where coefficients fi,Sj,rk depend on first two moments A and A 1 only, while all other 
coefficients w m , v n , u v depend just on first three moments A , A 1 and A 2 . 

A Hamiltonian structure of these integrable hydrodynamic chains is unknown at this 
moment. Just three Hamiltonian subcases have been found and completely investigated. 
The Hamiltonian hydrodynamic chains (here h 2y k = dh 2 /dA k , k = 0, 1, 2) 

A k = 2h 2i2 A k+l + h 2il A k + (k + 2)A k+ \h 2y2 ) x + (k + l)A k {h 2>1 ) x + kA k ~\h 2fl ) x 

are associated with the Kupershmidt-Manin Poisson bracket (see the second part in [9], 
[13] and [19]) 

{A k {x), A n {x')} = [kA k+n ~ l D x + nD x A k+n - x ]5{x - x'), k,n = 0, 1, 2, 

the Hamiltonian hydrodynamic chains (here = dh\/dA k , k = 0, 1) 

A k = (a + P)h hl A k+1 + l3h lfi A k + [a(k + 1) + 2f3]A k+ \h hl ) x + (ak + 2f3)A k (h lfi ) x 

are associated with the Kupershmidt Poisson brackets (see [8] and [17]) 

{A k {x), A n {x')} = [{ak + (3)A k+n D x + {an + (3)D x A k+n ]5{x -x'), k, n = 0, 1, 2, 

the Hamiltonian hydrodynamic chains (see [18]) 

A\ = (3ti {A°)A k+1 + {ak + 2P)A k+1 h^{A°)A° x 

are associated with the simplest Dorfman Poisson bracket (see [5]) 

{A k {x), A n {x')} = [{ak + (3)A k+n+1 D x + {an + (5)D x A k+n+1 ]5{x - x'), k, n = 0, 1, 2, ... 

If these hydrodynamic chains are integrable, then the corresponding Hamiltonian den- 
sities h {A°), hi(A°, A 1 ) and h 2 {A°, A 1 , A 2 ) cannot be arbitrary. The simplest case h (A°) 
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is described in [18] (see also [29]), while a full list of admissible expressions h^A^^A 1 ) 
and h2(A°, A 1 , A 2 ) is given in [8] and [9], respectively. 

Mainly, this paper is devoted to a description of integrable Hamiltonian hydrodynamic 
chains 

A* = 2(/M 1 + 6A° + OK(A°)Al + K(A°)(PAI + 5A% (3) 
A\ = /^°)(/3^ +1 + ^) + [(a^ k = 1,2,..., 

associated with more general Dorfman Poisson brackets (see [5]) 

{A k (x), A n (x')} = [T kn (A)D x + D x T n \A)]5(x -x'), k,n = 0, 1, 2, ... (4) 

where (a, (3, 7, S, e, £ are arbitrary constants) 

r 00 = PA 1 + 6A° + f , T kn = {ak + (5)A k+n+1 + (jk + 5)A k+n + ekA k+n -\ k + n > 0. 

Hamiltonian densities ho(A°) determining integrable hydrodynamic chains (3) are ex- 
tracted by virtue of the moment decomposition approach (see [29]). The crucial feature 
of this moment decomposition approach is a reducibility (i.e. the so called infinite dimen- 
sional analogue of the Darboux theorem) of Dorfman Poisson brackets to the canonical 
form d/dx. 

This paper is organized in the following way. In Section 2, the first (and most general, 
i.e. all constants a, /3, 7, 5, e, £ are arbitrary) class of integrable hydrodynamic chains as- 
sociated with Dorfman Poisson brackets is described. Corresponding Vlasov type kinetic 
equations are derived. Generating functions of conservation laws and commuting flows are 
found. A hierarchy of such integrable hydrodynamic chains is extended on negative mo- 
ments and negative time variables. Three dimensional two component hydrodynamic type 
systems (as well as three dimensional quasilinear equations of the second order) connected 
with this hierarchy are presented. In Section 3, the second class (i.e. /3 = but 5 7^ 0) of 
integrable hydrodynamic chains is described. In comparison with the general case, these 
hydrodynamic chains possess a momentum density, but Hamiltonian densities depend 
on infinitely many moments A k . In Section 4, the third (and the last, i.e. /3 = and 
5 = 0) class of integrable hydrodynamic chains is described. In Conclusion, a relationship 
between integrable Hamiltonian hydrodynamic chains and integrable Hamiltonian three 
dimensional two component hydrodynamic type systems is discussed. 

2 Dorfman Poisson brackets and the moment decom- 
position approach 

The above infinite component Dorfman Poisson brackets can be reduced to the canonical 
form (see [6]) 

{a i (x) : a k (x')} = — S'(x-x'), i,k = 1,2, ...,N, (5) 

where are arbitrary constants (8^ is a Kronecker symbol), due to the so-called "moment 
decomposition" 
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with an appropriate choice of functions fk.i{ al ) for any natural integer N. 

Theorem: Suppose ^ 0, in such a general case, Dorfman Poisson bracket (4) 
reduces to (5) under the moment decomposition 

dA k = (3- k J2 £ mV(a m )(V\a m )-S) k da m , fe = 0, 1, (7) 

where the function V(p) satisfies the ODE 

VV" = ^v"+(j- 2 -f)v' + Pe- 1 6 + 'f, (8) 

and £ is an integration constant of the simplest constraint 

±J2e m V 2 (a m ) = PA 1 + 6A + C (9) 

Proof: A substitution of moment decomposition (6) into Poisson bracket (4) implies 
the recursive relationships 

(ak + 0)fk+n+i,i + ilk + S)f' k+nji + ekf' k+n _ ltl = -f nti ftp k,n = 0,1,2,... 

and constraint (9). For k — 0, these ODE's 

PfUi,i + S fki = 7fn,if°,i> n = 0,1,2,... 
can be reduced to the common form ((3^0) 

£* = f?(| i -*) B > n = 0,l,2,... 

A substitution of these expressions into the above recursive relationships leads to (8), 
where V(a % ) = fhJe 1 . Moreover, taking into account (8), an integration of the differential 
d(aA k+2 + 7y4 fe+1 + eA k ) leads (see (7)) to an infinite series of constraints 

(ak + 2f3)A k+1 + ( 1 k + 25)A k + ekA k - 1 + t k = /T fc ^ e m \Z 2 (a m )(y'(a m ) -<5) fc , k = 0,1,..., 

where £ fc are integration constants. It is easy to see, if k — 0, then £ = 2£ (see (9)). The 
Theorem is proved. 

For any positive integer M, an arbitrary Hamiltonian density hM{A°, A 1 , A M ) and 
Poisson bracket (5) determine the Hamiltonian hydrodynamic type system (see [31]) 

i 1 a dh M 

a t = -d x — , (10) 
reducible to the symmetric form (see [27]; here hM, m = dh/dA m , m = 0,1, M) 

a t = ~ d * [ L^ 1 ""^" J = 9 * I V i a )l_, L h M,m I ■ (11) 

* \m=0 / \ m=0 / 



It is easy to verify that (11) is a hydro dynamic reduction of the Hamiltonian hydrodynamic 
chain (associated with Dorfman Poisson bracket (4)) 

M 

A° t = 2t{h Mfl ) x + Y}{un + 2P)A n+1 + { 1 n + 25)A n + enA n -\h M , n ) x 

n=0 

M 

+^[(cm + P) A x +l + (T« + S)K + enA n x - l ]h M ,n, 

n=0 
M 

a* = Yl^ k + n ) + w) Ak+n+l + di k + n ) + ^)A k+n + e(k + fii^-^y, 

n=0 
M 

+J2 [{an + + (771 + + enA k x +n - l ]h M>n 

n=0 

by virtue of moment decomposition (7). If the Hamiltonian density Hm {A , A 1 , A M ) 
is an arbitrary function, the above hydrodynamic chain is non-integrable, as well as its 
hydrodynamic reduction (11). Just in some special cases, this hydrodynamic chain and 
hydrodynamic reduction (11) became integrable. It means that the Hamiltonian density 
fiM{A°, A 1 , A ) must satisfy some overdetermined system in partial derivatives, which 
can be obtained utilizing different criteria of an integrability. For instance, in such a case, 
Haantjes tensor vanishes (see detail in [9]), a family of hydrodynamic reductions (see detail 
in [7]) parameterizes by iV arbitrary functions of a single variable. In this paper, integrable 
hydrodynamic chains are extracted by the method of symmetric hydrodynamic reductions 
(see detail in [27], [28], [29]). Without loss of generality, just the simplest Hamiltonian 
density ho(A°) is considered below. Indeed, if the Hamiltonian density h Q (A°) determines 
integrable hydrodynamic chain (3), then (3) possesses an infinite series of conservation 
law densities h n (A°, A 1 , A n ); then an infinite series of commuting flows (determined by 
these Hamiltonian densities h n (A°, A 1 , A n )) 

n 

A° tn+1 = 2£{h nfi ) x + K am + 2/3) A m+1 + ( 7 m + 25) A m + emA m ~ 1 } {h n , m ) x 

n 

+J2l( am + + iym + 5) A™ + emA™- 1 ]}!^, 

m=0 

(12) 

n 

A k n+1 = ^[(a(fc + m ) + ^)A k+m+1 + {y(k + m) + 25)A k+m + e(k + m)A k+m ^}(h ntm ) x 

m=0 

n 

+ J2i( am + P)A k x +m+1 + (im + 5)A k+m + emA k x +m - l ]h n , m 

m=0 

exists. Since M — 0, corresponding hydrodynamic reductions (11) simplify to 

4 = d x {v{d)h' Q ). (13) 
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The integrability criterion contains three steps only (see detail in [27], [28], [29] and 
[22]). 

1. Instead (13), let us introduce the generating function of conservation laws for 
hydrodynamic type system (13) replacing a* by p(x,t; A), i.e. 

Pti = d x {V{p)h' ), (14) 

where A is a free parameter (thus, a 1 = p(x,t, Aj), where A, are N distinct values of the 
parameter A). 

2. Under the semi-hodograph transformation p(x, t;X) —> X(x,t;p), (14) reduces to 
the Vlasov type kinetic (the collisionless Boltzmann type) equation (see [22]) 

X^ = V'(p)h' X x -V(p)X p h^A° x . (15) 

3. Finally, one can check (see [27]) the consistency (15) and corresponding hydrody- 
namic reduction (13). 

Since, all moments A k are expressed via field variables a n (see (7)), the last step can 
be replaced by a verification of consistency (3) and (15). In such a case, all necessary 
computations (see [8] and [9]) significantly simplify. Indeed, a substitution of X(p, A) into 
(15) leads to the linear PDE system (here X p = dX/dp, X k = dX/dA k , k = 0, 1, ...) 

°° ( h! \ 

VX P + J2i( ak + 2/3)^ fc+1 + (lk + 25)A k + ekA k ~ l ]X k + 2£ + (5 - V')-$ A = 0, (16) 



k=0 



A fc = ^, A; = 0,1,2,..., (17) 

where 

^ V -^- (18) 

A general solution of (17) is given by 

X(q,A) = B(q)J2— + C(q), (19) 

where B(q) and C(q) are not yet determined functions ("integration constants"). A 
substitution (19) into (16) leads to an extraction of an integrable hydrodynamic chain 
determined by the Hamiltonian density h (A°) such that (er is an arbitrary constant; 
a ^ 2/3) 



while the coefficients B(q) and C(q) in an asymptotic expansion (q — > oo) are given by 
their derivatives 

(Ingy- (a - 2(3)q 2 - 25q - e ^ ^ (aq-2^)B 



q(aq 2 + 7g + e) aq 2 + ^q + e 
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In a particular case a = 2/3, an integrable hydrodynamic chain is determined by the 
Hamiltonian density h (A°) such that (a is an arbitrary constant) 

ti = exp(A°/a), 

while the coefficients B(q) and C(q) in an asymptotic expansion (q — > oo) are given by 
their derivatives 

(hi BY- 25q + € C'= ^ aq ~ 2 ^ B 
1 ' q(2(3q 2 + iq + e) ' 2/3g 2 + 7 g + e' 

Equation (8) can be integrated in the parametric form (see (18)) 

J aq A + 7g + e J aq z + 7<? + e 

Then all conservation law densities hk can be found by a substitution of the inverse 
function g(A,A) (expanded in the Biirmann-Lagrange series, see, for instance, [20]) in 
(21) at the vicinity q —> oo. 

2.1 Negative conservation laws and commuting flows 

All positive commuting flows (12) possess infinite sets of conservation laws (cf. (1)) 

dtnhk = d x f nt k(h , hi, hk+n), k = 0, 1, 2, n = 2, 3, ... 

In the general case, integrable hydrodynamic chain (3) possesses local conservation laws 
for positive integers only. It means that any negative conservation law density must 
be expressed via infinitely many negative moments A~ n . Nevertheless, if e = 0, each 
negative conservation law density can be expressed via a finite number of negative 
moments A~ n . However, e = is not a particular case. Indeed, let us introduce new 
parameters Ck such that a = /?C2,7 — c i + 2<5c2,/3e = Co + c\5 + C2<5 2 . Thus, instead five 
arbitrary parameters a, (3,^,5, e, we shall use another five parameters 0, 5, Co, ci, c 2 . In 
such (8) simplifies to (see [28]) 

VV" = c 2 V' 2 + Cl V + c . (22) 

It means that integrable hydrodynamic chain (3) is parameterized by these three essen- 
tial parameters only, because the parameter 5 does not exist in generating function of 
conservation laws (14) as well as in the Hamiltonian density (see (20)) 

+ (23) 

while the parameter /3 can be incorporated to the integration constant a. It is easy to see 
below, that the parameter (3 can be removed from Poisson bracket (4) 

T kn = (3(c 2 k + \)A k+n+l + (jfe(ci + 25c 2 ) + 5)A k+n + °° + ° l6 + ^ kA k+n ~ l 
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by an appropriate scaling of moments A n . Thus, aforementioned integrable hydro dynamic 
chains are determined by one-parametric family of Dorfman Poisson brackets. Since 5 is 
a free parameter, let us choose the special value 5 of this parameter such that (i.e. e = 0) 

c + Cl 8 + c 2 t = 0. (24) 
Then integrable hydrodynamic chain (3) 

A° tl = 2(PA l + 5A + O(K) x + h' (PA l x + 5A x ), 

A% = [P(c 2 k + 2)A k+1 + (k( Cl + 28c 2 ) + 25)A%ti ) x + ti ((3A* +1 + 5A^ 
possesses an infinite series of negative local conservation laws 

d t ih- k = d x f k (h , h-!, h- k ), k = 0, 1, 2, ... 
The first negative conservation law is given by 

a^h^A- 1 ) = 5d x [h^(A^)K(A% (25) 

where 



, a \ 2S(1-C2)- C1 

h-i = ( A + 



26(1 - c 2 ) - ci 
Thus, the first negative commuting flow 

AU = h'_ i m-c 2 )A k x + (8(l-2c 2 )-c 1 )A k x - 1 } 

+ [f3(c 2 (k - 1) + 2)A k + (k( Cl + 25c 2 ) + 25(1 - c 2 ) - c 1 )A k ~ 1 } {h'_ 1 ) x 

is determined by the above Hamiltonian density /i_i(A~ 1 ). Then the "zeroth" conserva- 
tion law is 

dt-ihoiA ) = 13(1 - c^h'^A-^hoiA )). (26) 

Under moment decomposition (7) extended to negative values of index k, the above 
hydrodynamic chain transforms to the hydrodynamic reduction 



commuting with hydrodynamic type system (13). A corresponding generating function 
of conservation laws is given by (cf. (14)) 

"'—^(v^5 h '-^ 1) )- (27) 

Thus, obviously, generating function of conservation laws (14) is valid for Hamiltonian 
densities depended on negative moments too. However, a generating function of negative 



S 



conservation law densities is associated with another expansion (cf. (19)) at the vicinity 
q — > 0. Indeed, a similar computation as in the previous subsection leads to 

-l 

X(q, A) = B(q) + <?(?). (28) 



fc=— oc 

where 



n R V 2/3(c 2 -l)g + Cl + 2^(c 2 -l) ~, afi 
(m#J = -— ~= , C - 



/3c 2 g 2 + (ci + 25c 2 )g /3c 2 g 2 + (ci + 25c 2 )g 

Thus, a generating function of conservation laws for the integrable hydrodynamic chain 
determined by the Hamiltonian density hMi m 2 (A~ M2 , A' 1 , A , A Ml ) is given by (cf. 
(11), (14), (27)) 



(V'(p) - 6) 



m 



Pt = d x { V(p) ^ ^ h M t ,M 2 ,m 

\ m=-M 2 

where h MlM ^ m = dh MlM2 /dA m , m = 0, ±1, ±2, ... 

Remark: The hierarchy of integrable hydrodynamic chains determined by the Hamil- 
tonian densities h_ n (A~ n , A~ l ) and hk{A°, A k ) transforms to itself due to A n ~ x <H> 
A~ n and t n -H- t~ n . Indeed, let us introduce an auxiliary function V(p) such that (see 
(27)) 

V( p ) = V ^ (29) 
V'(p)-6 v ; 

This function V(p) satisfies (cf. (22)) 

VV" = c 2 V' 2 +c 1 V' + c , 

where 

25c 2 + ci 25c 2 + ci 

C 2 = : , Ci = 2(Co — 1)^-, : Co 

5(2c 2 -l) + ci' 12 ; 5(2c 2 -l) + Cl 

c = (c 2 - l) 2 2fc2 — + c 2 (l - c 2 ). 
b(2c 2 - 1) + ci 

Then an inverse transformation is given by 

V(p) = [5(1 - 2c 2 ) - Cl ] V(P) (30) 

l/'(p) + c 2 - 1 

It means, generating function of conservation laws (14) under the transformation t 1 -H- 
t -1 , A <H- A" 1 reduces to generating function of conservation laws (27). Thus, taking into 
account (29) and (30), all higher generating functions of conservation laws 

= d x [V(p)J2 d)m h k , m) , fc = 0,l,2,... (31) 

\ m=0 ^ / 

transform to corresponding lower generating functions of conservation laws 

Pt-» = d x (v(p) ^2 ^—J-^h-n,^ , n = 1, 2, ... (32) 

and vice versa. 
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2.2 Generating function of commuting flows and conservation 
laws 

These generating functions (31) and (32) can be incorporated in the sole generating func- 
tion of conservation laws and commuting flows 

d T{0 p(\) = d x G(p(\),p(0), (33) 

where the generating function of conservation law densities p(X) is determined above, while 
an auxiliary function p(Q is obtained from p(X) replacing A by (. Thus, a substitution of 
an expansion p(Q at the vicinity q — > oo or q — > into r.h.s. leads to generating functions 
(31) and (32) with an appropriate expansion of the so called "vertex" operator d T ^ with 
respect to parameter (. 

Theorem: The function G(p(X),p(Q) is defined by the quadrature 

Art m <r\\ ( n< V (p(Q) R( P (Q) \, mj V(p(X))R(p(Q) 

dG(p(\), P (0) = (Q(p(0) - v>(p(Q) - v'(p(\)) J M ) + vi(p(0)-v>( P W) Ma 

where 

V"(p(C)) 

R(p(0) = y^y> ^(P(0) = (C2 - 1)R(P(())- (34) 
Proof: The compatibility condition d t i(d T ^p(X)) = d T ^(d t ip(X)) (see (14)) implies 
d T{0 \nh' = Q(p(())d x \nh' + R(p(())d x p(0, (35) 
dG(p(X),p(0) , ( „ V(p(Q)R(p(Q) dG( P (X),p(0) V( P (X))R(p(0) 



dp(x) v(p(c))-v'(p(x)y dp(() v'(p(())-v(p(x)) 

where functions Q{p{()) and R(p(()) are not yet determined. The compatibility condition 

d dG(p(\),p(Q) _ d dG(p(\),p(Q) 
dp{X) dp(C) dp{() dp{X) 

leads to (34). Theorem is proved. 

If C2 = 1, then Q{p(C)) is a removable constant due to admissible shift in the vertex 
operator d T ^)+coiastd x — > <9 T (o- Then (35) reduces to the generating function of the 
"zeroth" conservation laws for all (positive and negative) commuting flows 

d T(O lnti = R(p(0)d x p((). 

It is easy to see ho = — lnh' in agreement with (23) for c 2 = 1. 
If c 2 7^ 1, then (35) reduces to (see the second equation in (34)) 

dr(oi(hT-~ 1 } = UQ(p(Cm' o r- 1 }- (36) 
It is easy to see h satisfies (23) for c 2 ^ 1. 



10 



2.3 Integrable three dimensional quasilinear equations of the 
second order 

The method of hydro dynamic reductions (see detail in [7]) allows to extract integrable 
three dimensional quasilinear equations of the second order (see [4]). In this subsection 
we present a list of some new such equations associated with the hierarchy of commuting 
hydrodynamic chains described above. 

Compatibility conditions (pt k )t n = {pt n )t k lead to integrable three dimensional hydro- 
dynamic type systems. If k = 1 and n — — 1, then the three dimensional hydrodynamic 
type system (see the general case in [23]) 

U t -i = (1 - c 2 )vu x - UV X , V t i = Suv x + (ci + (2c 2 - l)S)vu x 

is determined by the so called dispersionless Lax pair (see (14) and (27)) 

Ptl =d x (V(p)u), Pt- 1= ^( y,iy^ )» ( 3T ) 

where u = h' Q (A°) and v = h'_ 1 (A~ 1 ). If c 2 ^ 1 and £2 ^ 1 this hydrodynamic type system 
possesses two local conservation laws (25) and (26) 

(w C2-1 )t-i = (1 - c 2 )(vu C2 -%, (^ 2 " 1 ) tl = 5(uv £2 - l ) x . (38) 

Then the integrable quasilinear three dimensional equation of the second order (see [4]) 

ZtH-i z xt i = {z x ) c 2-^ dz xt -i H z xx (39) 

can be obtain introducing the potential function z such that z x = u c ' 2 ~ l and z t -i = 
(1 — C2)vu C2 ~ 1 . A similar equation can be derived from the second conservation law by 
another potential function z such that z x = v c ' 2 ~ l and z t i = buv C2 ~ x . In the particular 
case C2 = 1, (39) replaces by 

z tH -x = e Zx [Sz xt -i + (ci + (2c 2 - 1)5)^-1^], 

while the first conservation law in (38) replaces by (lnu) t -i = —v x . A substitution (where 
the potential function ip is introduced for (37)) 



into (38) implies the integrable quasilinear three dimensional equation of the second order 

Tp t i t -i + (Ci + C 2 5) WVt-^xw = K 1 - C V 7771 \ t l J t- lt l J xt 1 + 777777 ^a;; 



1 . 



If k = 1 and n — 2, the dispersionless Lax pair (see (14) and (31) where M — 1) 

Pfi = d x (V(p)u) , Pt2 = d x (V(p)w + V(p)V'(p)s) (40) 
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determines the three dimensional hydrodynamic type system (c 2 7^ —1) written in the 
conservative form 



c 2 + 1 V u c 2 J x c 2 + 2 

where u = h' (A°),w = h 10 — Shi^s = hi^. Moreover, the compatibility condition 
(Pt 1 )* 2 = (Pt 2 )* 1 leads to the constraint s = u C2+1 , which precisely coincides with a re- 
lationship between the corresponding Hamiltonian densities ho(A°) an d hi(A°,A l ), i.e. 
— {h'o) C2+l - Then the quasilinear three dimensional equation of the second order 

-Zxt? = z x z tH i + H {Zx)^- 1 z xt i H —{z x ) c 2-^z xx , (42) 



c 2 + 1 \c 2 — 1 c 2 / c 2 + 1 

is associated with the same (see (39)) potential function z of the first conservation law in 
(41). In the particular case c 2 = 1, (42) replaces by 

z xt 2 = 2z t i t i + 2(z t i + de Zx )z xt i + c e 2Zx z xx , 

while the first conservation law in (41) replaces by 

v ; \2u 2 ) x 
A substitution (where the potential function ip is introduced for (40); see also (37)) 

^ rr,f , S f 

u = T77-T' w = T77-T ~ V Wx) 



into (41) implies the integrable quasilinear three dimensional equation of the second order 



xt 2 



^1 + (C2+1) -P— ^i t l+Co ' 



The most interesting and exceptional case is given by c 2 = — 1. In such a case, 
dispersionless Lax pair (40) is no longer correct. The dispersionless Lax pair 

J¥ = d x {V(p)u) , Pt 2 = d x (V(p)w + V(p)V'(p)s + V(p)F{V'(p))) , (43) 

where F(q) satisfies 

qF(q) + (~q 2 + ciq + c )F'(q) = q 2 , (44) 
determines the pair of conservation laws (cf. (41)) 

o N / w + ci In u + Ci \ r . . 

(M" 2 ) t i+2 =0, u t 2 = w t i + co[w(lnu - I)]*. (45) 



Then the quasilinear three dimensional equation of the second order 

Z xt 2 = Z x Z tH i - + Ci(2 x ) 1/2 j Z xt i - y ln^e • 
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follows from the second conservation law of (45), while z is a potential function of the 
first conservation law of (45). 

A comparison (43) with (31) means that the Hamiltonian density h\ depends on A , A 1 
(as usual in the general case) and linearly on all other positive moments A k (see (12)). 
A corresponding expression for the function F(q) can be found by a substitution of the 
Taylor series F(q) = q 3 (e^ + e^q + e^q 2 + ...) into (44). Then all coefficients are determined 
iteratively, i.e. 

1 _ ci (fc - 3)e fc „ 2 - (fc - l)ci6 fc _i _ 

e 3 — o ' 64 ~~ _ T~2> 6fc — 7 ' ^ — 0,0,... 

The next integrable three dimensional quasilinear equation of the second order follows 
from (14) and (36) 

d T {u C2 ~ x ) = d x {Q{p)u c ^ x ), dtip = d x (V(p)u). 

Introducing a potential function z such that (utilizing the first conservation law) z x = 
m C2 ~ 1 , z T = Q(p)z x , the second conservation law transforms to 

1 V(v(s)) c 2 i 

z x z Tt i - z T z xt i = — - r —{z x )~ x z xx + (z^^V '{p{s))(z x z TX - Z T Z XX ), 

c 2 - 1 p'{s) 

where p(s) is an inverse function to Q(p) and s = z T /z x . Introducing a potential function 
xjj such that (utilizing the second conservation law) i/j x = p,ip t i = V(p)u, the first conser- 
vation law transforms to another integrable three dimensional quasilinear equation of the 
second order 



The most complicated integrable three dimensional quasilinear equation of the sec- 
ond order associated with an integrable hierarchy of commuting hydrodynamic chains 
presented in this Section (see (31) and (32)) can be obtained utilizing (33). The compat- 
ibility condition of two copies of (33) 

T ip(A) = d x G(p(\),p l ), <9 r2 p(A) = d x G(p(X),p 2 ) 

leads to the pair of conservation laws 

d T , P 2 = d x G{ P 2 y), d^p 1 = d x G( P \p 2 ), 

where d T i = St-(C)> &t 2 — dr^iP 1 = p((),P 2 = pijf) an d f] is an auxiliary parameter such £. 
Let us introduce two potential functions ip 1 and i\) 2 such that ip x = p 1 , ip T i = G(p 1 ,p 2 ) and 
ip x = p 2 , ip 2 i = G(p 2 ,p 1 ), then the above integrable three dimensional hydrodynamic type 
system reduces to two equivalent three dimensional quasilinear equations of the second 
order 
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where p 1 = G{i/>li,i/%) and p 2 = G(^ 2 ,^)- 

Under the transformation of independent variables x -H- t, generating function of 
conservation laws (14) 

Px = d t (V(p)u) 

reduces to the form 

ft = (!) > 

where p = f(s), s = V(p) and us = q. In the particular case c = 0, (22) reduces to the 
form (i.e. C2 = c, C\ = 1) 

V(p)V"(p) = cV'\p) +V'(p) 
by an appropriate scaling of the independent variable p. Then the dispersionless Lax pair 

* = «./(£), * = */(£), (46) 
determines the three dimensional hydrodynamic type system 

where 

s c — 1 

Also this three dimensional hydrodynamic type system 

U C ~ 2 CL C ~ 2 

a t = — ~{au x -ua x ), u y = — -(au x -ua x ) (47) 

u c — a c u c — a c 

can be rewritten as the three dimensional quasilinear equation of the second order (see 
the general case in [4]) 

2 V 2c (Z t ) „ , V c (Z t )-l 



Z x t H TTTTTT Z x Z ni 

tj 



Jxy (v c (z t )-iy xt ^ cv(z t ^ x 



where a = Z x and u = Z x /V(Z t ). Introducing a potential function ip such that (see (46)) 
q = ip x ,u = tp x /V(ip t ),a = ip x /V(ip y ), (47) reduces to the semi-symmetric form (i.e. this 
equation is invariant with respect to y -f->- 1; see the general case in [4]) 



V{%)V c {^ xt -V^ t )V c ^ y )^ xy 



This equation can be written in the form 
Let us take two extra copies of this equation 
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where r is the "fourth time" variable. Eliminating second derivatives ip XT , i/j xV x/j xy , one 
can obtain the remarkable three dimensional equation of the second order (see [4]) 

v e W> t )-v e W> y ) v c (i> y )-v c (j> T ) v^-v^i 
3 The first degenerate level /3 = 

In Section 2, the Theorem was formulated for the general case 7^ 0. This Section is 
devoted to this degeneration (3 = 0. 

Theorem: Suppose (3 = 0, but 5 ^ 0, in such a case, Dorfman Poisson bracket (4) 
reduces to (5) under the moment decomposition 

dA k = 5j2^ m a m {W'{a m )) k da m , k = 0,1, (48) 

where the function W(p) satisfies the ODE 

5 P W" = aW p + 7 W* + e, (49) 

and £ is an integration constant of the simplest constraint 

^e m (at = M° + e (50) 
Proof: Under the substitution V = (3W + Sp, (8) reduces to the form 

PWW" + 5pW" = aW' 2 + ~fW + e. 

Thus, the degenerate case (3 = is associated with (49). A substitution of moment 
decomposition (6) into the Dorfman Poisson bracket (cf. (4)) 

{A\x), A n (x')} = [F kn (A)D x + D x r nk (A)]5(x -x'), k,n = 0, 1, 2, ... 

where (a, 7, 5, e, £ are arbitrary constants) 

r 00 = SA° + f , T kn = akA k+n+1 + ( 7 fc + 5)A k+n + ekA k+n ~ l , k + n>0, (51) 

implies the recursive relationships 

anW k+n+l + {yn + 5)W k+n + enW k + n _ x = W k W' n , k,n = 0, 1, 2, 

where W n (a l ) = fni/ e i- Following the general case (see (7), (8) and below), suppose that 
a solution of this system is given by W m = W (W') m , where W (p) and W(p) are not 
yet determined. If n = 0, then W (p) = 5p (up to an additive constant). A substitution 
of the above ansatz into the recursive relationships leads to (49). Moreover, taking into 
account (49), an integration of the differential d(aA k+2 + "fA k+1 + eA k ) leads (see (48)) 
to an infinite series of constraints 

akA k+l + ( 7 fc + 25) A k + ekA k ~ l + £ fe = 5 2 ^ e m {a m ) 2 {W'{a m )) k , k = 0,1, 
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where £ fc are integration constants. It is easy to see, if k = 0, then £ = 2£ (see (50)). 
The Theorem is proved. 

For any positive integer M, an arbitrary Hamiltonian density fiM{A°, A 1 , ...,A M ) and 
Dorfman Poisson bracket (51) determine a hydrodynamic chain (see the previous Section), 
whose Hamiltonian hydrodynamic reduction (10) is presented in the symmetric form (see 
(48) and [27]; here h M , m = dh/dA m , m = 0, 1, M) 




5dJa*J2(W'(a l )) m h M , m 



m=Q 

Then such an integrable hydrodynamic chain possesses the generating function of conser- 
vation laws (cf. (31)) 

Pt = Sd x (p^(W'(p)) m h Mtt A . (52) 

V m=0 / 

In this case, the Hamiltonian density satisfies some overdetermined system. Following the 
general approach presented in the previous Section, we would like restrict our consider- 
ation to the simplest case ho(A°) only. However, this is impossible. In comparison with 
full (/3 7^ 0) Dorfman Poisson bracket (4), this reduced (j3 = 0) Dorfman Poisson bracket 
leads to the momentum density P = A , i.e. an arbitrary Hamiltonian hydrodynamic 
chain possesses a conservation law of the momentum 



MM M 

i_ a n— 1 s: i, 

lM 



A° t = 2£h M ,o + aY,nh M . n A n+l + ^( 7 n + 25)h M , n A n + e^tiM^" 1 - 5h } 



n=0 n=0 n=0 



It looks like the simplest integrable case is determined by the next Hamiltonian density 
/ii(A°,y4 1 ). However, corresponding hydrodynamic chain 

A° t = (2^1,0 + M 2 + ( 7 + 28)A 1 + eA°]h ltl + 2<L4%, - Sh) , 



A\ = [ a (k + l)A k+2 + ( 1 (k + l) + 25)A k+1 + e(k + l)A k }(h hl ) x 

+ [akA k+1 + (jk + 25) A k + ekA k - l }(hi, ) x + h^aA^ 2 + ( 7 + 5)A k+1 + eA k ] + 5h lfi A k 

depends on the highest moment A k+2 . Moreover, any higher commuting flow deter- 
mined by the Hamiltonian density h n (A°, A n ) depends on the highest moment A k+n+1 . 
Corresponding "time" variable t n+1 changes from n — 1, while t° must be reserved apri- 
ori for x, because the "zeroth" conservation law density ho(A°) is a momentum density 
A , which cannot create an "intermediate" hydrodynamic chain containing the high- 
est moment A k+1 . The momentum density generates just "trivial" commuting flow, i.e. 
A k = A k . This is a reason to identify x and t°. Nevertheless, a commuting flow deter- 
mined by the "time" variable t 1 exists. 
Lemma: The hydrodynamic chain 

A% = [2(£ + 8A°)f" + 5f + eh}A° x + [(7 + S)h + 2eb 2 \A\, (53) 
A% = [akA k+1 +{ 1 k+25)A k +ekA k - 1 }f'Al^5f+eb 1 )A k x ^{^ k = 1,2, ... 



1(3 



is determined by the Hamiltonian density depended nonlinearly on the "zeroth" moment 
and linearly on all higher moments, i.e. 



K = f(A°) + Y^nA" 



(54) 



n=l 



where b\ and b 2 are arbitrary constants, while all other constants b n satisfy the recursive 
relationships 



a(n - l)& n _i + (jn + 6)b n + e(n + l)fo n+ i = 0, n = 2, 3, 



(55) 



Proof: can obtained by a straightforward calculation. 

Moreover, integrable hydrodynamic chain (53) is associated with the generating func- 
tion of conservation laws (see (52)) 



where 



In this case (see [22]), 



p t i=5d x (pf(A°) + W(p)), 

oo 

W(p)=pJ^b m (W'(p)) m . 

m=l 

p w" = c 2 w' 2 + aw' + c , 



(56) 
(57) 
(58) 



where c^ are some constants. 

Theorem: Integrable hydrodynamic chain (53) is determined by Hamiltonian density 
(54), where 

en, a o\ (7 + S)h + 2e6 2 



(59) 



+ tr[( 7 + 5)6! + 2eb 2 ] ' 

a is an integration constant, and constants bk satisfy (55). 

Proof: Under the semi-hodograph transformation p(X, x, t) — > X(p,x,t), (56) reduces 
to the Vlasov type kinetic equation (see [22]) 

V = S[(f(A°) + W'(p))X x -pX p (f%}. 

A substitution X(p, A) leads to 



SpX p 



(7 + S)h + 2eb 2 ^ «A;A fc+1 + ( 7 Jfe + 25) A k + e/U*" 1 



/"(A ) 

and (17), where (instead (18)) 



fc=0 



An 



(60) 



SW'(p) - eh 



(7 + 5)6! + 2e6 2 ' 
Taking into account (see (58) and (61)) 

5 2 pW"(p) 



(61) 



(7 + 5)h + 2eb 2 



c 2 q 2 + dq + c , 



(62) 
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a substitution (19) into (60) implies the constraints c 2 = a, c,\ — 7, c = e, (59) and 

q(aq 2 + 7<? + e) aq + 79 + e 

Taking into account (61) and the constraints C2 = a, ci = 7, Co = e, a comparison (49) 
with (62) implies q = W'(p), i.e. 



(7 + 5)6i + 2e6 2 



Differentiating (57) and eliminating W from the above relationship, finally, (56) reduces 
to 

W = dxWf + + ((7 + S)b. + 2eb 2 )W), (64) 

where all constants bk satisfy (55). Theorem is proved. 
Equation (49) can be integrated in the parametric form 

/Sd/Q f 
— — — , W = qdp (65) 
aq 2 + 75 + e J 

Then all conservation law densities hk can be found by a substitution of the inverse 
function g(A,A) (expanded in the Burmann-Lagrange series, see, for instance, [20]) in 
(65) at the vicinity q — > 00. 

3.1 Negative conservation laws and commuting flows 

If hydrodynamic chain (53) is integrable, then an infinite series of higher (positive) com- 
muting flows (55) exist whose corresponding Hamiltonian densities h n (A°, ...,A n ) depend 
on a finite set of moments only. 

If e 7^ 0, hydrodynamic chain (53) does not possess negative local conservation laws, 
i.e. any negative conservation law density h- n depends on all negative moments A~ k . 
However, (49) is reducible to the desirable form 

5pW" = aW' 2 + (2ac + j)W' 

due to the shift W = W + cp, where the shift constant c is a solution of the quadratic 
equation ac 2 + 7c + e = (cf. (24)). Then reduced Dorfman Poisson bracket (51) 
transforms to the more simple form (expressed via new moments A k ) 

{A k (x), A n (x')} = [r kn (A)D x + D X Y nk {A)]5(x -x'), k,n = 0, 1, 2, ... 

where 

f 00 = SA° + e, T kn = akA k+n+1 + [(7 + 2ac)k + 5}A k+n , k + n > 0, 

and A = A , A 1 = A 1 + cA°, A 2 = A 2 + IcA 1 + c 2 A°, A 3 = i 3 + 3cA 2 + 3c 2 A 1 + c 3 A° etc. 
This is a consequence following from comparison (see (48)) 

dA k = 5Yjtma m {W'{a m )) k da m = 5 ^ e m a m (H>'(a m ) + c) k da m 
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with 



dA k = 5^e m a m (W\a m )) k da r 



Thus, without loss of generality, the choice e = allows to simplify all computations 
in this Section. In such integrable hydro dynamic chain (53) 

A% = [2(£ + 5A°)f" + 5f']A° x + ( 7 + *)Mi, (66) 

A% = [akA k+l + ( 7 fc + 25)A k ]f"A° x + Sf'A k + ( 7 + 5)hA k x + \ k = ±1, ±2, ... 

possesses also an infinite series of lower (negative) local conservation laws. For instance, 
the first negative conservation law is given by 



where 



while (see (59)) 



d t ^ 1 (A- 1 ) = Sd x if'(A°)h_ 1 (A- 1 )], 

5 



(lnfc-i)' = 
f"(A° 



(25 - rfA' 1 + a(j + 5)h 
( 7 + S)bx 



aA° + cr( 7 + 8)bi 
Thus, the first negative commuting flow 

A k ., = [ a (A;-l)A fe +( 7 (A;-l)+2(5)A fe - 1 ](/il 1 ) a .+/i / „ 1 [-a^+(5- 7 )^.- 1 ], k = 0, ±1, ±2, ... 

is determined by the first negative Hamiltonian density h_i(A~ l ). Its first non-negative 
conservation law (see (54)) is 

d t -iK = d x [b 1 ( 1 + 5)A hi l - a A f'hi l -5f'h^ 1 + (25- 1 )f'A~ 1 h'_ l ] 



+d x 



h '-i ( 5 - 7)$>^ n_1 - «5>^ n 



n=l 



n=l 



where constants b n are given by the reduced relationships (see (55)) 

anb n + ( 7 (n + 1) + 5)b n+1 =0, n — 1, 2, ... 

Under moment decomposition (48) the above hydrodynamic chain transforms to the hy- 
dro dynamic reduction 

Thus, the generating function of conservation laws for the first negative commuting flow 

p '-'=^G4) ft '-') (67) 

is a particular case of (14). Indeed, if e = 0, then ordinary differential equations (22) and 
(49) are connected by the transformation 

V(jp)W'(p)=p, 
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where (see (21) and (65); b is an integration constant and b 7^ 0) 

V(p) = V - (ftp"* - a) , W(p) = 7 / , f , (68) 
7 v ' J bp s — a 

and (22) reduces to the particular case 

While a computation of positive conservation laws is based on expansion (19) at the 
vicinity q — > 00, negative conservation laws can be found utilizing expansion (28) at the 
vicinity q — > 0. A generating function of conservation law densities 

7 



V — [a + 



q(X,A) 



is an inverse expression to q = W'(jp) (see (68); without loss of generality, constant b can 
be fixed to the unity). The first series of higher (positive) conservation law densities can 
be found by a substitution of an inverse expansion g(A, A) from (19), where (see (63)) 

(lnBV= a<1 ~ 26 C '= a ^ + ^ biq ~ 2 ^ B q^oo 
q(aq + 7)' q(aq + 7) 

The second series of lower (negative) conservation law densities can be found by a substi- 
tution of an inverse expansion q(X, A) from (28), where 

q{aq + V {aq + l)q 

Then substitution of these series of conservation law densities to the generating function 
of conservation laws (64) 

Pti =d x {5fp+{ 1 + 5)b l W) (70) 

allows to extract both series of conservation laws. 

Remark: Functions B(q) and C(q) are found by the same computation as in the 
previous Section. The Vlasov type kinetic equation 



V = [Sf + (7 + 5)bxW]X x - 6p\ p f"A° x 



is connected with (70) by a semi-hodograph transformation p(x,t; X) — > X(x,t;p). Taking 
into account the negative part of integrable hydrodynamic chain (66) 

A% = [akA k+1 + ( 7 fc + 2S)A k }f"A° x + Sf'A k x + (7 + 5)hA k x + \ k = -1, -2, 

a substitution of expansion (28) into this Vlasov type kinetic equation leads to (69). 
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3.2 Generating function of commuting flows and conservation 
laws 

An integrable hierarchy of commuting hydro dynamic chains described in this Section can 
be embedded into a sole generating function of conservation laws and commuting flows 
(33). The function G(p(X),p(()) can be found from the compatibility condition of (33) 
and (70) 

d t i P (X) = d x [up{\) + W(p(X))}, d T{0P (X) = d x G(p(\),p(()). 
Theorem: The function G(p(X),p(Q) is defined by the quadrature 

where 

8W"(v(0) a 

R(p(0) = . Q(p(0) = f fl(p(0). (7i) 

Proof: The compatibility condition d t i(d T ^p(X)) = d T ^(d t ip(X)) implies 

d T (c)U = Q(p(())u x + R(p(0)d x p(0, (72) 
dG(p(X),p(0) -Q(p( C) ). P(C)fl(p(0) dG( P (X),p(0) P(X)R( P (0) 



d P (x) ^ vrvsy; w'(p(0)-w'ip(x)y d P (o w>(p(0)-w>(p(x)y 

where functions Q{p{()) and R(p(()) are not yet determined. The compatibility condition 

d dG(p(\),p(Q) d dG(p(\),p(Q) 
dp{X) dp(C) dp{() dp{X) 

leads to (71). Theorem is proved. 

Moreover (72) reduces to (see the second equation in (71)) 

d T(0 e au/S = (Q(p(0)e au/ % 

It is easy to see the conservation law density e au ^ s is nothing but a momentum density 
A (up to unessential additive and multiplicative constants). 

3.3 Integrable three dimensional quasilinear equations of the 
second order 

Compatibility conditions {ptk)t n = {Pt n )t k l ea d to integrable three dimensional hydrody- 
namic type systems. If k — 1 and n — — 1, then three dimensional hydrodynamic type 
system 

ct 7 — 5 

Ut- 1 + v t o + —vu t o = 0, v t i = uv t a H — vu t o (73) 

o o 

is determined by the dispersionless Lax pair (see (67) and (70)) 

P 

W'(p) 



Pti = d t o(up + W(p)), pt-i = d t o ( — ) , (74) 
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where 

u = I i jru /'» u = 
(7 + 

Three dimensional hydrodynamic type system (73) can be written in the conservative 
form (7 7^ 5) 



(e w/5 )r. + -j(ve au / s ) t0 = 0, (u^) t i = (uv^) t o. 



OL is 5 5 

Introducing the potential function z such that (see the first conservation law) 



5 Sz t -i 

u—— hi Zip, v — , 

a az t o 

(73) reduces to the three dimensional quasilinear equation of the second order 
a(z t oz t -i t i - z t -iz t o t i) = 5z t o \nz t o ■ z t -i t o + (7 - 5 - 5\nz t o)z t -iz t o t o; 
introducing the potential function z such that (see the second conservation law) 

V = (Ztp) * , u = — , 
Z t o 

(73) reduces to another three dimensional quasilinear equation of the second order 
5(z t az t i t -i - z t iz t o t -i) + (z t o)V[( 7 _ S)zto - az t i]z t o t a + a(5 t o) V +1 5 4 o 4l = 0. 
A substitution (where the potential function ip is introduced for (74)) 

u = , v = lp t -l 

into (73) implies the integrable quasilinear three dimensional equation of the second order 

(aP^(^o) + (7 - 5)^o - a^i) ^-i^o t o = 0. 



If k = 1 and n = 2, then the dispersionless Lax pair (see (52)) 

p t i = d t o (up + W(p)), p t 2 = d t o (wp + spW'(p)) (75) 
determines the constraint s = e au ^ s and the three dimensional hydrodynamic type system 

(e au/s ) t i = (w+(u- H*±I ] e Q "/ 5 j , u t2 + ™ ( o = + um t o, (76) 



V a / /to 

which is equivalent to the three dimensional quasilinear equation of the second order 

aSz t o t 2 = a 2 z t oz t i t i+a('y+S—25\iiZto)z t oz t o t i+S[5zto(ln z t o) 2 — (7+25) z t o (In ^0— l)+a^i]^oto, 
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where z is a potential function for the conservation law in (76). A substitution (where 
the potential function ip is introduced for (75)) 

fa-WiM ip t 2 , faap t i - aW(ip t o) 

u = ■ , w = W {ip t o) exp ' 



lp t o lp t o V 

into (76) implies the integrable quasilinear three dimensional equation of the second order 

exp I fafl ) ' r'"' 2 ~ ~ ^ tOta ) = a ^ tH1+ \ 7 + ^ ) ^ tH1 



I ^2 + 5 ) ^ <W (W ) W- 

The next integrable three dimensional quasilinear equation of the second order follows 
from (70) and (72) 



d T e au/s = d t o(Q(p)e au/s ), d t ip = d t o(up + W(p)). 

Introducing a potential function z such that (utilizing the first conservation law) z t o = 
e au/s^ z ^ _ Q(pjzip, the second conservation law transforms to 

$P( S ) i Tin/ i \ 

Z t oZ Tt i - Z T Z t o t i = —Z t oZ t o t o + {—\nZfi + W (p(s)))(z t oz Tt o - Z T Z t o t o), 

ap{s) a 

where p(s) is an inverse function to Q(p) and s = z T /z t o. Introducing a potential function 
ip such that (utilizing the second conservation law) ipt° = P^t 1 = up + W(p), the first 
conservation law transforms to another integrable three dimensional quasilinear equation 
of the second order 

4 The second degenerate level (3 = and 6 = 

In Section 3, the Theorem was formulated for the case /3 — but 6^0. This Section is 
devoted to the second degeneration 6 = 0. 

Theorem: Suppose (3 = and 6 = 0, in such a case, Dorfman Poisson bracket (4) 
reduces to (5) under the moment decomposition 

dA k = ^e rn {U\a m )) k da m , A; = 0,1,..., (77) 

where the function U (p) satisfies the ODE 

U" = all' 2 + 7C/' + e, (78) 
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and £ is an integration constant of the simplest constraint 

( 79 ) 

Proof: Under the substitution V = /3U + 5p + 1, (8) reduces to the form (cf. (49)) 

PUU" + (6p + l)U" = aU' 2 + jU' + e 

Thus, the double degenerate case (3 = and then 5 = is associated with (78). A 
substitution of moment decomposition (6) into the Dorfman Poisson bracket (cf. (4)) 

{A k (x),A n (x')} = [T kn (A)D x + D x T nk (A)}5(x-x'), k, n = 0, 1, 2, ... (80) 

where (a, 7, e, £ are arbitrary constants) 

r oo = ^ r fcn = + 7 ^fc+« + e ^+"-i) ; jfe + n > 0, 

implies the recursive relationships 

n(aU k+n+1 + 7£4 +n + eC/ fc+n _i) = U k U' n , k, n = 0, 1, 2, 

where t/ n (a*) = f' ni l^i- Following the general case (see (7), (8) and below), suppose that 
a solution of this system is given by U m = Uo(U') m , where Uo(p) and U(p) are not yet 
determined. If n = 0, then Uo{p) = 1 (an arbitrary nonzero additive constant is fixed to 
the unity here). A substitution of the above ansatz into the recursive relationships leads 
to (78). Moreover, taking into account (78), an integration of the differential d(aA k+2 + 
jA k+l + eA k ) leads (see (77)) to an infinite series of constraints 

k{aA k+1 + 1 A k + eA*- 1 ) +U = Y, ^m{U\a m )) k , k = 0, 1, 

where £ fc are integration constants. It is easy to see, if k = 0, then £ = £ (see (79)). The 
Theorem is proved. 

For any positive integer M, an arbitrary Hamiltonian density Iim{A , A M ) and 
Dorfman Poisson bracket (80) determine a hydrodynamic chain (see the previous Sec- 
tions), whose Hamiltonian hydrodynamic reduction (10) is presented in the symmetric 
form (see (77) and [27]; here h M , m = dh/dA m ,m = 0, 1, M) 

* \m=0 / \ni=0 

Then such an integrable hydrodynamic chain possesses the generating function of conser- 
vation laws (cf. (31)) 

Pt = d x (^(t/'b)) m V m ) • (81) 

\m=0 / 

In this case, the Hamiltonian density satisfies some overdetermined system. Since this 
double degenerate case ((3 = and 5 = 0) is very familiar to the degenerate case (,5 = 0), 
very similar results are presented below. 
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Lemma: The hydrodynamic chain 

A% = {2if + eM° + (761 + 2eb 2 )A 1 ) x , (82) 

A k tl = k(aA k+1 + 1 A k + eA k ~ l )f"A x + eb x A\ + (761 + 2eb 2 )A k+ \ k = 1,2, ... 

is determined by Hamiltonian density (54) where b\ and b 2 are arbitrary constants, while 
all other constants b n satisfy the recursive relationships 

a(n — l)& n -i + l~nb n + e(n + l)b n+1 = 0, n — 2, 3, ... (83) 

Proof: can obtained by a straightforward calculation. 

Moreover, integrable hydrodynamic chain (82) is associated with the generating func- 
tion of conservation laws (see (81)) 

Pti=d x (f\A°) + U(p)), (84) 

where 

00 

Uip) = J2 b m(U'(p)) m . (85) 

771=1 

In this case (see [28]), 

U" = c 2 U' 2 + Cl U' + co, (86) 

where are some constants. 

Theorem: Integrable hydrodynamic chain (82) is determined by Hamiltonian density 
(54), where 

f{A) ~ a + aAO ' (87) 

a an integration constant and bk satisfy (83). 

Proof: Under the semi-hodograph transformation p (x, t; A) — > X(x,t;p), (84) reduces 
to the Vlasov type kinetic equation (see [22]) 

A tl = U'(p)\ x - \ p f"(A°)A x . 

A substitution X(p, A) leads to 

x '=[ nA°) " 2f -^ ? j A » (88) 



and (17), where (instead (18) and (61)) 



U'(p) - eh 

q = — — -i. 89 

y -fh + 2eb 2 y J 



Taking into account (see (86) and (89)) 

U" 



761 + 2eb 2 
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c 2 q +Cig + c , (90) 



a substitution (19) into (88) implies the constraints c 2 = a, c.\ — 7, c = e, (87) and 

(lnB)>= ( f~ e v C= a 2 q ~ 2i B. (91) 
q{aq 2 + 7g + e) aq l + 7g + e 

Taking into account (89) and the constraints C2 = a, c\ = 7, Co = e, a comparison (78) 
with (90) implies q = U'(p), i.e. 



761 + 2eb 2 ' 



Differentiating (85) and eliminating U' from the above relationship, finally, (84) reduces 
to 

p t i = d x {f'{A°) + (761 + 2eb 2 )U(p) + eb lP ), (92) 

where constants bk satisfy (83). Theorem is proved. 

Equation (78) can be integrated in the parametric form 



dq U= f 

aq A + 75 + e ' J aq 2 + 75 + e 



P = / r..., ' / 2 ,y ', y • (93) 



Then all conservation law densities hk can be found by a substitution of the inverse 
function g(A, A) (expanded in the Biirmann-Lagrange series, see, for instance, [20]) in 
(93) at the vicinity q — > 00. 

4.1 Negative conservation laws and commuting flows 

If hydrodynamic chain (82) is integrable, then an infinite series of higher (positive) com- 
muting flows exist whose corresponding Hamiltonian densities h n (A°, A n ) depend on a 
finite set of moments only. 

If e 7^ 0, hydrodynamic chain (82) does not possess negative local conservation laws, 
i.e. any negative conservation law density /i_ n depends on all negative moments A~ k . 
However, (78) is reducible to the desirable form 

U" = all' 2 + (7 + 2ac)U' 

due to the shift U = U + cp, where the shift constant c is a solution of the quadratic 
equation ac 2 + 7c + e = (cf. (24)). Then reduced Dorfman Poisson bracket (80) 
transforms to the more simple form (expressed via new moments A k ) 

{A k (x), A n (x')} = [T kn (A)D x + D x f nk {A)]6(x -x'), k,n = 0, 1, 2, ... 

where 

f 00 = f , Y kn = k(aA k+n+1 + (7 + 2ac)A k+n ), k + n > 0, 

and A = A , A 1 = A 1 + cA°, A 2 = A 2 + 2cA l + c 2 A°, A 3 = A 3 + 3cA 2 + 3c 2 A 1 + c 3 A° etc. 
This is a consequence following from comparison (see (77)) 

dA k = J2 e m (U'(a m )) k da m = Y^e m {U'{a m ) + cfda m 
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with 

dA k = J2 e m(U'(a m )) k da m . 

Thus, without loss of generality, the choice e = allows to simplify all computations 
in this Section. In such integrable hydro dynamic chain (82) 

A% = (2f /' + lhA 1 ) x , A% = k(aA k+1 + iA k )f"A° x + 7 M* +1 , k = ±1, ±2, ... (94) 

possesses also an infinite series of lower (negative) local conservation laws. For instance, 
the first negative conservation law is given by (7 7^ 0) 

d t i ln( 7 A" 1 - a) = -—d x ln(aA° + a) (95) 

a 

where (see (87)) 

f\A Q ) = ^MaA Q + a). 
a 

Thus, the first negative commuting flow 

AU = (k- l)(aA k + 1 A k ~ l )(h'_ l ) x - h!_ x {aA k x + 1 A k ~ l ) 

is determined by the first negative Hamiltonian density h_i(A~ 1 ) = In^A^ 1 — a). Its 
first non- negative conservation law (see (54)) is 

(2 00 
-t-h ln(a + aA°) + 7 M°/i'_i - (aA° + 1 A- 1 )fh'_ 1 - Y)> n (aA n + -yA^ti^ 

where constants b n are given by the reduced relationships (see (83)) 

anb n + 7 (n + 1)6„ +1 = 0, n = l,2, ... 

Under moment decomposition (77) the above hydrodynamic chain transforms to the hy- 
drodynamic reduction 

Thus, the generating function of conservation laws for the first negative commuting flow 

»-' = 9 -m- (96) 

is a particular case of (14). Indeed, if e = and 5 = 0, then ordinary differential equations 
(22) and (78) are connected by the transformation 

V(p)U'(p) = 1, 

where (see (21) and (93)) 

V(p) = , U(p) = - In 7 - - ln(l - ae 7P ) (97) 

7 a a 
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and (22) reduces to the particular case 

VV" = V (V + a) . 

While a computation of positive conservation laws is based on expansion (19) at the 
vicinity q — > oo, negative conservation laws can be found utilizing expansion (28) at the 
vicinity q — > 0. A generating function of conservation law densities 

p — in I a + 



7 V ?(A,A) 

is an inverse expression to g = t/'Qo) (see (97)). The first series of higher (positive) 
conservation law densities can be found by a substitution of an inverse expansion g(A, A) 
from (19) (see (91)) 

A(g,A) = (xg-2^1ng + (ag + 7)2j — , q oo. 

k=0 1 

The second series of lower (negative) conservation law densities can be found by a substi- 
tution of an inverse expansion q(X, A) from (28) 

-l 

A(g,A) = (ag 2 + 7 g)^g- fc - 1 ^-(Tg, q 0. (98) 

fc=— oo 

Then substitution of these series of conservation law densities into generating function of 
conservation laws (92) 

761 aA° + a 
= d x In — 99 

a 1 — ae^ p 

allows to extract both series of conservation laws (the above expression is obtained by 
virtue of (97)). 

Remark: Expansion (98) is found by the same computation as in the previous Sec- 
tions. The Vlasov type kinetic equation 



A, = t*. (U'MK - A p -^) 



is connected with (99) by a semi-hodograph transformation p(x,t; A) — > X(x,t;p). Taking 
into account the negative part of integrable hydrodynamic chain (94) 

A% = k{aA k+1 + jA k )f"A° x + jhA k x +1 , k = -1, -2, 

a substitution of expansion (28) into this Vlasov type kinetic equation leads to (98). 
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4.2 Generating function of commuting flows and conservation 
laws 

An integrable hierarchy of commuting hydro dynamic chains described in this Section can 
be embedded into a sole generating function of conservation laws and commuting flows 
(33). The function G(p(X),p(()) can be found from the compatibility condition of (33) 
and (92) 

dtipiX) = d x [u + U(p(\))}, d T(oP (X) = d x G(p(X),p(()). 
Theorem: The function G (p(\) , p(Q) ^ s defined by the quadrature 

*W), P (0) = (q«0) - u^XiM) ) * (A » + m SU) W ' 

where 

R(p(0) = ex V [2aU(p(()) + 7P(C)], Q'(p(0) = <*R(p(Q) (100) 
Proof: The compatibility condition d t i(d T ^)p(X)) = d T ^(d t ip(X)) implies 

d T{0 u = Q(p(())u x + R(p(0)d xP (0, (101) 
dG(p(X),p(0) - Q{p{C) y R (P(0) dG(p(X), P (0) R( P (0) 



dp(x) u'( P (0)-u'(p(x)Y dp(c) u'(p(())-u'(p(x)) 

where functions Q(p(C)) and R(p(()) are not yet determined. The compatibility condition 

d dG(p(\),p(Q) d dG(p(\),p(Q) 
dp(X) dp(C) dp{() dp{X) 

leads to (100). Theorem is proved. 

Moreover (101) reduces to (see the second equation in (100)) 

d T{0 e au = (Q(p(())e au ) x 

It is easy to see the conservation law density e au is nothing but a momentum density A 
(up to unessential additive and multiplicative constants). 

4.3 Integrable three dimensional quasilinear equations of the 
second order 

Compatibility conditions {ptk)t n = {pt™)t k l ea d to integrable three dimensional hydrody- 
namic type systems. If k = 1 and n — — 1, then three dimensional hydrodynamic type 
system 

ik- 1 + v t° + Oivu t o =0, v t i = jvuto (102) 
is determined by the dispersionless Lax pair (see (96) and (92)) 

p t i = d t o{u + U(p)), p t -i = d t o (rj^j , (103) 
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where 

1 7 

u — — ln(aA° + cr N 



a 1 — a 

Hydro dynamic type system (102) can be written in the conservative form 

(Ot-i + a(t;e Q > = 0, Qnv) t i = W . 

Introducing the potential function z such that (see the first conservation law) 

1 , z t -i 
u = — In z t o , v — , 

a az t o 

(102) reduces to the three dimensional quasilinear equation of the second order 

a(z t oz t -i t i - z t -iz t o t i) = jz t -iz t o t o; 
introducing the potential function z such that (see the second conservation law) 

S n Z t l 

v = e ■ *°, u = — , 
7 

(102) reduces to another three dimensional quasilinear equation of the second order 

z tH -i + e 5 ' (nfZfifi + aztp t i) = 0. (104) 
A substitution (where the potential function ip is introduced for (103), see (97)) 

u = $ + _ ln(l - ae^% v = 1 ^ 
a 1 — ae w *° 

into (102) implies the integrable quasilinear three dimensional equation of the second 
order 

+ i ^T^t-i(a^to t i + 7^*0*0 ) = 0. 

1 — ae 7V «° 

If /c = 1 and n = 2, then the dispersionless Lax pair (see (81) and (99)) 

Pti = d t o(u + U(p)), Pt 2 = d t o(w + sU'{p)) (105) 

determines the constraint s = e au and the three dimensional hydrodynamic type system 
written in the conservative form 

a(e au ) t i = (aw - ^e au ) t o : u t 2 = w t x. (106) 

Introducing the potential function z such that (see the first conservation law) 

1, 7 

u — — In Zto , w = z t i H — z t o, 
a a 

(106) reduces to the three dimensional quasilinear equation of the second order 

z t o t 2 = z t a(az tH i + 72^1 ); 
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introducing the potential function z such that (see the second conservation law) 

U = Z t i, W = Z t 2, 

(106) reduces to another three dimensional quasilinear equation of the second order 

z t o t 2 = e^* 1 (az tH i + jzfi t i). (107) 
A substitution (where the potential function ip is introduced for (105), see (97)) 

u = 4> t i + -ln(l - ae 7 ^' ), w = ^ 2 - ie ^ +a ^ 
a 

into (106) implies the integrable quasilinear three dimensional equation of the second 
order 

Vv^ = e^ 1 (1 - «e^)( 7 ^o tl + a^). 

Remark: The whole hierarchy of positive and negative commuting flows possesses the 
symmetry t k+1 -H- t~ k . Without loss of generality it is enough to prove for the independent 
variables t 1 -h- t° only Indeed, first generating function of conservation laws (92) 

Pt = d x (u + U(p)) 

can be rewritten in the same form 

q x = d t (v + U(q)), 

where q = U(p),p = U(q),v t = —u x (see (95)). Since (see (78) where e = and (97)) 
U"(p) = aU'\p) + 7*7'(p), then U"(q) = - 7 t/' s '' (q) - aU'(q). Thus, equations (104) and 

(107) are equivalent to each other (by an appropriate change of the constants a 7 and 
the independent variables t 1 <H- t°,t 2 <H- 1^ 1 ). 

Let us consider the three dimensional hydrodynamic type system written in the con- 
servative form (see (99), (101) and (97)) 

ie au / 1 

d T e au = d t o 1 d tl p = d t0 (u--\n(l-ae^) 

1 — ae^ p \ a 

(108) reduces to the integrable three dimensional quasilinear equation of the second order 

a(z t oz Tt i - z T z t o t i) = (z T - 7^0) z Tt o, 

where the function z is a potential of the first conservation law; (108) reduces to another 
integrable three dimensional quasilinear equation of the second order 

VVti - 7^*o 4 i = e 7 ^°(77/vto + a^ Tt i), 
where the function if; is a potential of the second conservation law. 
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5 Conclusion 



In this paper, integrable Hamiltonian hydrodynamic chains associated with Dorfman Pois- 
son brackets are described. These Dorfman Poisson brackets are parameterized by six 
(at, (3, 7, 5, e, £) constants. However, three distinguish cases are selected by special values 
of two constants (3 and 5 only. The first and most general class is determined by arbitrary 
parameters (a- ^ 0, (3 ^ 0, 7 7^ 0, 5 7^ 0); the second class is determined by the restriction 
(3 = while 5^0; the third class is determined by more deep restriction (3 = and 5 = 0. 
In the particular case 5 = 0, hydrodynamic chains of the first class were constructed 
in [27]; corresponding Hamiltonian reductions are the Hamiltonian chromatography sys- 
tems. Hydrodynamic chains of the second class also are well known. In a particular case 
(a special limit a = 0), they are nothing but the Kupershmidt hydrodynamic chains (see 
[3], [17], [21] and [26]). The most investigated Benney hydrodynamic chain (see [1], [12], 
[13], [14], [15], [16], [19], [32]) belongs to the third class. 

In a general case, hydrodynamic chains of the third class are considered in [28]. How- 
ever, the simplest hydrodynamic chain (82) was missed. This hydrodynamic chain is found 
(as well as corresponding hydrodynamic chain (53) of the second class) in the presented 
paper. Corresponding Hamiltonian density depends on infinitely many field variables 
(moments A k ). This phenomenon was never mentioned in a literature. 

All integrable hydrodynamic chains possess infinite series of conservation laws and 
commuting flows. Corresponding generating functions are constructed. These hydro- 
dynamic chains are extended on negative values of moments A k as well as associated 
integrable hierarchies are extended on negative values of time variables t n . 

Most obvious and simple three dimensional hydrodynamic type systems are presented. 
They are converted to the three dimensional quasilinear equations of the second order (see 
a general theory in [4]). Most of them are new. Nevertheless, the list of these equations 
(as well as associated hydrodynamic type systems) is so large, that all already found 
examples in [10] can be extracted from this list by virtue of different parametric reductions. 
Moreover, the list presented in [10] can be split on three parts according given approach 
here. Examples (18) and (21) belonging to the Benney hierarchy are of the third class; 
examples (20)i, (20)2 and (22) belong to the Kupershmidt hierarchy (see [26], where 
(3 = 3); most complicated examples (16), (17) and (19) belong to the first class, where 
V(p) = p(p) and the elliptic Weiershtrass function satisfies p' 2 (p) = 4p 3 (p) + 5 2 , while 
the function V(p) satisfies 2V(p)V"(p) = 3(V'\p) - 5 2 ) (see (22)). This is a full list of 
integrable three dimensional Hamiltonian hydrodynamic type systems (see detail in [10]) 




In another paper [11], a full list of integrable three dimensional Hamiltonian hydrodynamic 
type systems 

( u \ ( d x \ f h u \ 
\vj t \d. d y )\h v ) 

also is presented. Examples (21) and (24) belong again to the Benney hierarchy; examples 
(20), (22)x and (23) belong to the Kupershmidt hierarchy (see [26], where (3 = 3); example 
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(22) 2 belongs to the first class, where V(p) = pip). And just the last example (19) is 
not yet recognized. Thus, most these examples associated with two dimensional constant 
Poisson brackets simultaneously are connected with hydrodynamic chains associated with 
Dorfman Poisson brackets (but not vice versa!). However, the question "how to connect 
both Hamiltonian structures" is open at this moment. 

The relationship between three dimensional hydrodynamic type systems, three di- 
mensional quasilinear equations of the second order and hydrodynamic chains described 
above is very important. All hydrodynamic chains (2) presented in this paper are very 
convenient for more deep investigation by the moment decomposition approach (see de- 
tail in [29] and in [24]). This approach allows to extract multi-parametric solutions (see 
detail in [30]) of these hydrodynamic chains and corresponding three dimensional hydro- 
dynamic type systems as well as three dimensional quasilinear equations of the second 
order. Thus, any aforementioned example can be equipped by a corresponding hydrody- 
namic chain described in this paper. It means, that such a hydrodynamic chain can be 
considered as an infinite set of the so called pseudo-nonlocalities (i.e. moments) allowing 
to extend integrable three dimensional two component hydrodynamic type systems on 
infinitely many field variables. It means that a complexity of integrable three dimensional 
quasilinear equations of the second order can be translated to a complexity of two dimen- 
sional hydrodynamic chains (i.e. hydrodynamic type systems containing infinitely many 
equations) . 
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